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Abstract. In the Scheduling Machines with Capacity Constraints prob- 
lem, we are given k identical machines, each of which can process at most 
mi jobs. A'l jobs are also given, where job j has a non-negative process- 
ing time length tj > 0. The task is to find a schedule such that the 
makespan is minimized and the capacity constraints are met. In this 
paper, we present a 3-approximation algorithm using an extension of 
Iterative Rounding Method introduced by Jain [4]. To the best of the 
authors' knowledge, this is the first attempt to apply Iterative Rounding 
Method to scheduling problem with capacity constraints. 
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1 Introduction 

We consider the Scheduling Machines with Capacity Constraints problem (SMCC): 
There are k identical machines, and machine i can process at most nii jobs. Given 
M < X]i<i<fc jobs with their processing time lengths, we are to find a sched- 
ule of jobs to machines that minimizes the makespan and meets the capacity 
constraints. 

Scheduling problem is a classical NP-Hard problem and has been stud- 
ied extensively. In the general setting, we are given set T of tasks, number 
k of machines, length l{t,i) G Z+ for each t ^ T and machine i £ 
the task is to find a schedule for T, namely, a function f : T ^ [^--k], to 
minimize maxjgji jt] ^^^^ ^^^^^^ Z(t, i) . Lenstra, Shmoys and Tardos [6] gave 
a 2- approximation algorithm for the general version and proved that for any 
e > no (§ — e)-approximation algorithm exists unless P = NP . Their method 
based on applying rounding techniques on fractional solution to linear program- 
ming relaxation. Gairing, Monien and Woclaw [2] gave a faster combinatorial 
2-approximation algorithm for the general problem. They replaced the classical 
technique of solving the LP-relaxation and rounding afterwards by a completely 
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integral approach. For the variation in which the number of processors k is 
constant, Angel, Bampis and Kononov [1] gave a fully polynomial-time approxi- 
mation scheme (FPTAS). For the uniform variation where /(t, i) is independent 
of the processor i, Hochbaum and Shmoys [3] gave a polynomial-time approxi- 
mation scheme (PTAS). 

The SMCC problem is one of the uniform variations, with capacity con- 
straints on machines. One special case of SMCC problem in which there are 
only two identical machines was studied in [8] [10] [11]. Woeginger [9] gave a 
FPTAS for the same problem. General SMCC problem is a natural generaliza- 
tion of scheduling problem without capacity constraints and can be used in some 
applications in real world, such as students distributions in university, the Crew 
Scheduling problem in Airlines Scheduling [12] [13], etc. In the Crew Scheduling 
problem, crew rotations, sequences of flights legs to be flown by a single crew 
over a period of a few days, are given. Crews are paid by the amount of flying 
hours, which is determined by the scheduled rotations. Airline company wants 
to equalize the salaries of crews, i.e. to make the highest salary paid to crews 
minimum. Rotations starts and ends at the same crew base and must satisfy a 
large set of work rules based labor contracts covering crew personnel. In the con- 
cern of safety issues, one common contract requirement is the maximum times 
of flying of a single crew in a period of time. So the aim is to find a scheduling 
of rotations to crews that minimizes the highest salary and meets the maximum 
flying times constraints. 

In many literature, researchers approached scheduling problem using round- 
ing techniques. Lenstra, Shmoys and Tardos [6] applied rounding method to the 
decision problem to derive a p-relaxed decision procedure and then used a binary 
search to obtain an approximation solution. In the SMCC problem, the capacity 
constraints defeat many previous methods. In this paper, our algorithm is one of 
the class of rounding algorithms, but use a different rounding method introduced 
by Jain [4]. We do not round off the whole fractional solution in a single stage. 
Instead, we round it off iteratively. 

Iterative Rounding Method, introduced by Jain [4], was used in his break- 
through work on the Survivable Network Design problem. This rounding method 
does not need the half-integrality, but only requires that at each iteration there 
exist some variables with bounded values. In [4], Jain observed that at each it- 
eration one can always find a edge e has Xe at least 1/2, which ensures that the 
algorithm has an approximation ratio of 2 . As a successful extension of Jain's 
method, Mohit Singh and Lap Chi Lau [7] considered the Minimum Bounded 
Degree Spanning Trees problem and gave an algorithm that produces a solution, 
which has at most the cost of optimal solution while violating vertices degrees 
constraints by 1 at most. As far as the authors know. Iterative Rounding Method 
has been used in graph problems, and has produced many beautiful results. 

In this paper, we apply Iterative Rounding Method to the scheduling problem 
with capacity constraints and obtain a 3-approximation algorithm. To the best of 
the authors' knowledge, this is the first attempt to approach scheduling problem 
with capacity constraints using Iterative Rounding Method. 
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The rest of the paper is organized as follows. In Section 2, we formulate the 
SMCC problem as an Integer Program, give its natural relaxation and introduce 
our relaxation, Bounded Linear Programming Relaxation (BLPR). In Section 
3, we present some properties of BLPR and prove theorems that support our 
algorithm. In Section 4, we present bounding theorems and an approximation 
algorithm, IRA, and prove that it has an approximation ratio of 3 . 

2 Preliminary 

Formally, the SMCC problem is as follows: Given a positive integer k, k positive 
integers {mi\mi > 0,1 < i < k}, M non-negative integers {tj\tj > 0,1 < j < 
M < '^i}' ^6 to solve the following Integer Program (IP): 

l<i<k 

l<i<k 
l<j<M 
k,l<j<M 

There are some relaxations, one of which is the following natural Linear 
Programming Relaxation (LPR) dropping the integrality constraints. 

l<i<k 
l<i<k 
l<j<M 
k,l<j<M 

We don't use LPR directly, but use an alternative relaxation. Bounded Linear 
Programming Relaxation (BLPR): Given a positive integer k, k positive integers 
{mi\mi > 0,1 < i < fc}, M non-negative integers {tj\tj > 0,1 < j < M < 

X]i=i "^i}) ^ vector b = (&i, 62, • • • , ^fe) and C {(i,j)\l < i < k,l < j < 
M}, find a feasible solution under the following constraints 

l<i<k (1) 

l<i<k (2) 

l<j<M (3) 

ihj)^^ (4) 
l<i<k,l<j <M (5) 



mmimize c 
subject to '^jLi Xijtj — c < 

EM ^ 
j—l Xij \ JTTj 



G {0,1} 1 < i < 



minimize c 
subject to ^ij^j — c<0 

EM ^ 
j—l Xij 2i l^i 



Xij >0 1 <i < 



M 



subject to Xijtj < bi 

M 

^ ^ Xij ^ TTlj 



k 



1=1 



Xij — 1 

Xij > 
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where vector 6 = (61, 62, • • • , &fe), called upper bounding vector, is added to de- 
pict the different upper bounds of machines more precisely, and is added to 
represent the partial solution in algorithm. Each £ !F indicates that job 

j has been scheduled to machine i . Those {xij\{i,j) G J^} are considered as 
constants. 

We will show that properly constructing vector b = (61, &2, • • • , bk) makes the 
solution produced by our algorithm under control and easy to analyze. 

Definition 1. In a BLPR problem A, upper bounding vector b = (61, 62, ... , bk) 

is called feasible if A is feasible. 

Keeping the upper bounding vector b feasible all the time is the key of our 
algorithm, which guarantees that we can always find a feasible solution bounded 
by 6 . 



3 Techniques 



Before we present our algorithm, we need to introduce some properties of BLPR. 

With respect to the partial solution T, let Ci denote the number of already 
scheduled jobs in machine i, namely, q = G ^}\ ■ Note that rui — Ci 

indicates the free capacity in machine i . 

We call a job free if it has not been scheduled to any machine and call a 
machine free if it still has free capacity. For a feasible fractional solution x to 
A, define a bipartite graph G(x) = G{L, R, E), called supporting graph, where 
L represents the set of free machines, R represents the set of free jobs and 
E = {{i,j)\xij > 0,{i,j) ^ T} . We denote the number of free jobs and the 
number of free machines by M* and k* respectively. Note that for free job j, 

Consider the Constraint Matrix of A, which consists of the coefficients of the 
left side of equalities and inequalities, except for the non-negativity constraints 
from (5): 



/(tit2 ... ^m) 



(11... 1) 



(ti t2 ■■ ■ Im) 



(11... 1) 



1 1 



V V 



V V 



(11... 1) 
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1 
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where the 1®^ to A:* rows represent the constraints from (1), the {k + 1)^ to 
(2A;)*^ rows represent the constraints from (2), and the (2fc + l)*^ to (2fc + M)*^ 
rows represent the constraints from (3). 

One can verify that the (fc + 1)*^ row can be hnearly expressed by the rest of 
rows. Thus the rank of Constraints Matrix is bounded by the following lemma 

Lemma 1. Constraints Matrix has a rank at most M + 2k — 1 . □ 

Recall that, a basic solution x to Tl is the unique solution determined by a 
set of linearly independent tight constraints that are satisfied as equalities. We 
remove all zero variables in x so that no tight constraints comes from (5). Thus 
the number of non-zero variables in x never exceeds the rank of Constraints 
Matrix. When =$, the following inequality holds 



We can remove those non- free machines from A, move fixed variables {xij G 
to the right side of the equalities and inequalities as constants and remove 
variables fixed to 0. By doing this, we obtain a new sub-problem and only focus 
on free jobs and free machines. In the new sub-problem. Lemma 1 holds. So in 
general, the following corollary holds. 

Corollary 1. Given a BLPR problem, A, its basic solution x and supporting 
graph G{x) = G{L,R,E), we have 



We introduce lemmas on the basic solution to A when there are no less free 
jobs than twice the free machines, namely, M* > 2k*. 

Lemma 2. // A is feasible with M* > 2k* and x is a basic solution, there exist 
M* — 2k* + 1 variables with values of 1 . 

Proof. For a basic solution x, we construct supporting graph G(x) = G{L, i?, E) . 
Suppose that I of M* free jobs are of degree of at most one in G. Note that each 
of them has degree at least one. Each of the rest M* — I free jobs has degree of 
more than one. The following inequality holds. 



< M -I- 2fc - 1 



(7) 



\E\ < M* +2/c* - 1 



□ 



\E\ >2{M* -l) + l = 2M* -I 



(8) 



By Corollary 1, we have 



1>M* -2k* + 1 



(9) 



□ 



The following corollary holds 



Corollary 2. // A is feasible with M* > 2k* and x is a basic solution, there 
exist a free machine p and a free job q such that Xpg = 1 . □ 
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4 A 3-approximation Algorithm 

In this section, we present an approximation algorithm IRA . Let A denote the 
makespan in the solution produced by IRA, OVT denote the makcspan in the 
optimal solution to A . 

We introduce three bounding theorems on BLPR. Noting that when M* > 
2k* we can find a Xpq = 1, we can schedule g to p without increasing the lengths 
in the fractional solution. We first show the theorem for the case M* > 2k*. 

Theorem 1. Given a BLPR problem, A, with M* > 2k* and its basic solution 
X . Based on A, we construct a new BLPR problem f{A) as follows: 

1. Find a variable Xpq = 1 ; 

2. r^J'+{p,q) ; 

3. The rest parts of f{A) are the same as A . 

If b is a feasible upper bounding vector of A then b' is a feasible upper 
bounding vector of f{A) . □ 

When M* < 2k* and some free machines have free capacity of one, we have 

the following theorem. 

Theorem 2. Given a BLPR problem, A, with M* < 2k* and its basic solu- 
tion X. Moreover some free machines have free capacity of 1. Based on A, we 
construct a new BLPR problem g{A) as follows: 

1. Let p denote a m,achine with 1 free capacity; 

2. Let q denote the free job with the largest length; 

3. r^J^+ip,q) ; 

I 6' ^(61,62,..., + tq,bp+i, ...,bk) ; 

5. The rest parts of g{A) are the same as A . 

If b is a feasible upper bounding vector of A then b' is a feasible upper 
bounding vector 0/5 (^) . 

Proof. To schedule q to p, for each free machine s ^ p with x^q > 0, we move 
Xgq fraction of job q to machine p then move back from p as much as possible 
but no more than x^q fraction of jobs other than q as in Algorithm 1. 

Because q has the largest length among the free jobs, we can guarantee that 
the length of each free machine p' ^ p will not increase and the length of machine 
p will increase by at most tq . Note that in gi^A) , Xpq = 1 and p is no longer 
free. This implies there is a feasible solution to gi^A) . □ 

When M* < 2k* and every free machine has more than 1 free capacity, we 
can schedule these M* jobs arbitrarily but only assuring that each machine gets 
no more than 2 jobs. One can prove the following theorem. 

Theorem 3. Given a feasible BLPR problem,. A, with M* < 2k*. Moreover 
every free machine has free capacity of at least 2. We construct a new BLPR 
problem h{A) as follows: 
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Algorithm 1 G — transition 

Require: A BLPR A with M* < 2k* and machine p has free capacity of 1 
Ensure: A BLPR g(A) 

1: Let X be a basic feasible solution to A and q be the longest free job; 

2: while there exists p' p such that Xp/g > do 

3: if there exists q' q such that Xpg/ > then 

4: a^m.m{Xp/g,Xpg>} ; 

6- -^pq' * '^pq' ) 

7: else {Vg' # q, Xpg, = 0} 

8. ^pq' ^ ^pq' ^p'q ) 

9: end if 
10: end while 



1. Schedule free jobs arbitrarily hut only assuring that each machine gets no 
more than 2 jobs; 

2. For each machine p, increase bp by the sum of job (s) scheduled to p and 

update T accordingly; 

3. The rest parts of h{A) are the same as A . 

By doing this, we obtain h{A) in which all jobs have been scheduled and h{A) is 
feasible . □ 

By Theorem 1, as long as M* > 2k*, we can always schedule a free job to a 
free machine but without increasing its length. If M* < 2k*, Theorem 2 and 3 
guarantee wc still can make our decision in a fairly simple way. We present our 
algorithm using Iterative Rounding Method, IRA, in Algorithm 2. 

Finding a basic solution to a linear program can be done in polynomial time 
by using the ellipsoid algorithm [5] then converting the solution found into a 
basic one [4] . Together with the following observation 

Lemma 3. At Line 3, y is a feasible \ipper bounding vector of A . □ 
the correctness of IRA follows from Theorem 1, 2 and 3. 

Corollary 3. Algorithm IRA always terminates in polynomial time. □ 

The analysis of the performance of IRA is simple with the help of upper 
bounding vector b, noting that once a component of b is increased, the machine 
will be no longer free. We now show that IRA is a 3-approximation algorithm. 

Theorem 4. IRA is a 3-approximation algorithm. 

Proof. Consider any machine p with length A in the solution produced by IRA . 

Note that once bp is increased at Line 9 machine p will no longer be free. So 
exactly one of the following statements is true when the algorithm terminates: 

1. bp hasn't been increased, then bp < yp ; 

2. bp has been increased once at Line 9, then bp < yp + tq for some q ; 
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Algorithm 2 IRA 
Require: An IP A 

Ensure: A feasible integral solution 
1: Construct natural linear programming relaxation F ; 

2: Solve r optimally and let y = (j/i, J/2, ■ ■ ■ , 2/fc) be the lengths of machines in the 
optimal solution; 

3: Construct a BLPR A, letting y be the upper bounding vector and = ; 
4: while M* > do 
5: if M* > 2k* then 

6: A^f{A) ■ 
7: else {M* < 2k*} 

8: if there exists a machine p with free capacity of 1 then 

9: A^g(A) ■ 

10: else {every free machine have more than 1 capacity} 

11: A <- h{A) ■ 

12: end if 

13: end if 

14: end while 

15: return ; 



3. bp has been increased once at Line 11, then bp < yp + tq^ + tg^ for some 

Note that yp and maxg {tq} are two trivial lower bounds of OVT . Also note 
that after the algorithm terminates, the integral solution produced by IRA, 
contained in J^, is also bounded by upper bounding vector b . By definition of 
feasible upper bounding vector, we have inequality 

A<bp< 30VT (10) 

as expected. □ 

5 Conclusion 

In this paper, we consider the SMCC problem, a uniform variation of general 
scheduling problem, which has capacity constraints on identical machines. Using 
an extension of Iterative Rounding Method introduced by Jain [4], we obtain a 
3-approximation algorithm. This is the first attempt to use Iterative Rounding 
Method in scheduling problem and it shows the power of Iterative Rounding 
Method. It is still luiknown that whether the approximation ratio can be im- 
proved or whether the Iterative Rounding Method can be used to obtain a good 
approximation algorithm for the non- uniform version of scheduling problem with 
capacity constraints. 
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